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Abstract: In the event of an observation of neutrinoless double beta decay, a relevant 
question would be: what lepton number violating physics is responsible for the decay? 
The exchange of Major ana neutrinos and/or supersymmetric particles may contribute. We 
point out that measurements of supersymmetric signals at the LHC, including single slepton 
production, could be used to help bound some supersymmetric processes contributing to 
neutrinoless double beta decay. LHC information about the supersymmetric spectrum 
could be combined with B^-B^ mixing data in order to bound a competing neutrinoless 
double beta decay process involving sbottom exchange. 
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1. Introduction 



Neutrinoless double beta decay {0^(3(3) provides the most sensitive probe of lepton number 
violation as well as a Majorana nature of neutrinos. At the quark level, ^vj3j3 corresponds 
to the simultaneous transition of two down quarks into two up-quarks and two electrons. 
While the most prominent mechanism in the literature that causes this decay involves the 
exchange of a massive Majorana neutrino, several other possibilities have been discussed. 
Here, we focus on the attractive alternative where Qv(3(3 is mediated by the exchange of 
sparticles in supersymmetric models with R-parity violation [1-5]. 

R-parity violating couplings arise in a general supersymmetric extension of the Stan- 
dard Model (SM), where the superpotential contains renormalisable baryon- and lepton- 
number violating operators. The presence of both sets of operators typically leads to vi- 
olation of stringent bounds on proton decay [6], unless the parameters are unnaturally 
suppressed. Proton decay bounds are evaded if a discrete symmetry is imposed which 
forbids at least one set of such parameters. A widely-studied example is R-parity [7], un- 
der which both sets of operators are odd under the parity transformation and hence are 
forbidden. This symmetry also has the advantage of having the lightest supersymmetric 
particle (LSP) as a natural dark matter (DM) candidate. On the other hand, there ex- 
ist R-parity conserving (RPC) dimension 5 operators which could potentially lead to fast 
proton decay [8,9]. One way to suppress these operators is by instead imposing proton 
hexality [11,12]. 

Instead of R-parity, we focus on an alternative, namely baryon triality (also known 
as baryon parity in some literature) [10, 13]. This Z3 discrete symmetry allows for the 
dimension 4 R-parity violating (RPV) terms which violate lepton number, while those that 
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violate baryon number are forbidden. A significant advantage of this class of models is that 
dimension 5 proton decay operators are forbidden. The LSP will decay via the non-zero 
RPV couplings present, hence a neutralino LSP cannot be a dark matter candidate. Other 
dark matter candidates are viable, for example the gravitino, since its decay is slow on 
cosmological time scales [14] . Alternatively the dark matter could originate from a hidden 
sector. In the rest of this paper, we use the term R-parity violation to denote the lepton 
number violating, R-parity violating interactions. A survey of effective lepton number 
violating operators may be found for example in [15]. 

After imposing a symmetry that forbids baryon-number violating terms, the RPV 
superpotential is 

yVRPV = 2^ijkLiLjEk + X'^jj^LiQjDk — HiLiHu, (1.1) 

where we have suppressed all gauge indices and used the notation of Ref. [16] . The {i, j, k} G 
{1, 2, 3} are family indices. Ajj^, \'-^ are dimensionless trilinear RPV couplings, and jJLi are 
bi-linear RPV parameters, having dimensions of mass. 

Compared with the RPC minimal supersymmetric extension to the Standard Model 
(RPC MSSM), the presence of additional RPV couplings leads to distinctive signatures 
at a collider and have interesting physical consequences, for example providing neutrino 
masses. As the RPV couplings violate lepton number by 1 and involve lepton doublets, 
they automatically lead to Majorana neutrino masses [17-25] without the need to introduce 
additional field content such as a right handed neutrino. Qv(3f3 proceeds through the '1-1' 
entry in the neutrino Majorana mass matrix mpp which connects two electron neutrinos in 
a basis where the charged lepton mass matrix is diagonal. We refer to this mechanism as 
ruf^p contribution in the remainder of the paper. For three left-handed Majorana neutrino 
masses mj, with PMNS mixing matrix Uai, where a G {e,//, r} and i G {1,2,3}, 

3 

"^/3/3 = XI "^i^er (1-2) 

mpi3 may be a complex quantity. There are additional higher dimensional effective opera- 
tors, characteristic of this class of models, which mediate {)v(3(3 without an mpp insertion. 
As these operators violate lepton number by two units without the need of an mpp inser- 
tion, they are not suppressed by the smallness of the neutrino masses as a result. These 
channels will be called direct contributions. In general, the direct and rrij^p contributions 
contribute to Ovf3f3 simultaneously, so depending on their relative magnitudes, interference 
between the direct and nipp decay matrix elements may need to be included in determining 
the decay rate of Quf3f3. 

A measurement of the Oz^/?/? rate alone does not fix a neutrino mass scale, since it is 
possible that direct contributions are non-negligible. However, the RPV couplings leading 
to the direct contributions can affect other observables, which could then be used to con- 
strain the amplitude of the direct contributions. If they can be experimentally bounded, 
one may infer mpfj from the 0i'f3/3 decay rate. 
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The aim of this paper is to explore the interplay between direct and mp^ contributions, 
in particular how different constraints and observations may shed light on the underlying 
mechanisms of ^v(3(3. We will show how the combined knowledge of the masses of the 
electron sneutrino and sbottoms and constraints from B^-B^ mixing could allow us to 
determine an upper bound on a direct Qv(3(3 channel involving sbottoms. Searches for 
single slepton production at the LHC could provide valuable information on the value of 
which (with measurements of various sparticle masses) will allow one to bound its 
direct contribution to Qiy(3(5. A first exploration on the relationship between OvPP decay 
rate and single slepton production at the LHC may be found in [26]. A complementary 
way to probe mechanisms of Oiy(3f3 was discussed in [27, 28] by combining measurements 
from different nuclei. Efforts to relate neutrino masses and collider phenomenology in other 
theories with lepton number violation may be found for instance in [29]. 

For concreteness, we only consider Qvf3f3 of ^^Ge. In the rest of this paper, numerical 
values of nuclear matrix elements (NMEs) and half life all refer to this nucleus. A related 
work on contributions of trilinear RPV terms to Qvf3f3 is presented in [30]. There, particular 
attention is paid to nuclear matrix element calculations and contributions from different 
sparticles in the presence of a non-zero A'^^^^^ coupling. We go beyond the scope of this 
work in several ways: most importantly, we focus on what may be inferred from different 
experimental measurements on the mechanism that produces Oi^/3/3. We have also corrected 
certain terms in the effective Lagrangian at the quark level, that were incorrect in the 
literature. 

This paper is organised as follows. In section ^ we briefly review the possible OvPP 
mechanisms in the RPV MSSM. The current experimental half life limit Ty^^ of ''^Ge, as 
well as the neutrino oscillation data are summarised in section ^. We also list some useful 
scaling relations between the RPV parameter bounds and SUSY breaking parameters used 
there. In section 4.1 we proceed to discuss how the RPV contribution to B^-B^ mixing can 
affect the possibility of the A']^^3A'^3^ direct contribution to be the dominant observable dvf5f3 
channel. We then investigate the prospects of observing the single selectron resonance at 
the LHC and its implication of Qvf3[3 in section 4^, before concluding in section Technical 
information about parton-level OuPP calculations is in Appendix^and our NMEs are listed 
in Appendix ^ 



2. Mechanisms of Oiy(3p in the RPV MSSM 

In the RPV MSSM, it is possible to construct Majorana neutrino mass models that explain 
the neutrino oscillation data [13,21,25,31-36]. Oi^/3/5 could proceed through standard light 
Majorana neutrino exchange with the m/3/3 mass insertion. In addition, there are direct 
contributions via the RPV coupling products A'^]^3A']^3]^ and A'j^^j^A'^]^;^. 

The RPV coupling products A^^^A^-.^^ contribute to the neutrino mass matrix {m^)ij, in 
particular they generate the '1-1' entry m^^. A Feynman diagram using the mass insertion 
approximation (MIA) [37] with ^n^^i^i is shown in fig. [l|. In this approximation, m^^ can 
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be written as 



■»..^|^ y"'1" ln(%+(^»^). (2.1) 

OTT^ mi —mi \mi J 



Stt^ ml — ml \ml 

oll ojifi bjifi 

Here m? , m? and m? represent the entries of the sbottom mass matrix in an ob- 
vious notation, and m^ is the running mass of the down quark. Ov^f] may then proceed 
via exchange of a virtual Majorana neutrino with a insertion. A Feynman diagram 
depicting this process is displayed in fig. ||. For a realistic model, one expects there to be 
many non- vanishing bi- linear and/or tri-linear RPV operators present in order to fill out 
the effective 3x3 neutrino mass matrix with non-zero entries. Bi-linear RPV couplings 
could also lead to direct Oi^/3/5. However, as discussed in [39], neutrino mass terms obtained 
from these couplings that are consistent with the observed neutrino mass scales typically 
lead to negligible direct contributions to the Oz^/?/3 rate. On this basis we neglect their 
direct contribution, but bear in mind the possibility that they may enhance beyond 
what is expected from the tri-linear couplings. We also expect that the Xijk couplings may 
contribute to m^^, but not to affect direct {)uf3f3 significantly. Coupling products of the 
form y^'iik^'iki and X'^^nXnk could violate lepton number i and electron number by 1 unit 
each. For i 7^ 1, their contributions to direct via PMNS mixing are suppressed by 

the mass scale of the light neutrinos. As a result their contributions to direct Oi^/3/3 is likely 
to be subdominant and will not be discussed further in this paper. 

{rnf)Ln 
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Figure 1: An example diagram showing a contribution to mpj3 from the product A']^]^3A']^3]^. 
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Figure 2: A Feynman diagram showing 0i/f3/3 via exchange of a virtual Majorana neutrino. 
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A Feynman diagram representing a direct contribution mediated by X'n^Xi^i is shown 
in fig. 1^. As the corresponding matrix element does not contain neutrino mass insertions, 
it is not suppressed by the smahness of the neutrino mass scale. In principle, the product 
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Figure 3: Direct A']^]^3A']^3]^ contribution to Qv(3(3. 



'^'ii2'^'i2i could also lead to Oi^pp. However, this coupling product is tightly constrained by 
K^-K^ mixing [38], and hence this contribution is neglected in the rest of this paper. 

It should be noted that A'^^^^A'^^^^^ can also mediate Oz//3/3 via a diagram similar to fig. |3|. 
However, in the case of A'j^^j^A'^j^^, there are other diagrams contributing to OvP/S which 
dominate. These diagrams are shown in fig. 0. 
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Figure 4: A'j^j^j^ A'j^j^j^ contributions to Qvj3l3. 
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3. Experimental limits 



The most stringent lower limit on the Ge Oi//3/3 half life was bounded by the Heidelberg- 
Moscow experiment [40,41] to be 

T°/f ^ > 1.9 • lO^Syrs. (3.1) 

The most common interpretation of such a limit is in terms of a model in which only mpis 
contributes to Qv(3(3. In the event of an observation in the next round of OvPP searches, 
the half-life will be used to infer |m/3/3|, as in fig. ^. The figure shows that in the event 
of an observation of Ty^^ < 10^^ years, Impf^] > 50 meV. |m/3/3| > 10 meV implies a 
non-hierarchical (i.e. inverted or quasi-degenerate) spectrum [42]. 
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Figure 5: Value of mpp inferred from a future observation of 
only Majorana neutrino exchange contributes. 
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for our NMEs and assuming 



Such inferences depend crucially on the assumption that no other process (for example, 
direct RPV processes) contribute to 01^(3(3. We shall now take into account the possibility 
that the direct RPV processes can simultaneously contribute to m^^. Depending on the 
RPV couplings considered, different mechanisms dominate the Oi//3/3 process. Thus eq. ( |3.1| ) 
can be translated to upper bounds on particular products of RPV couplings. For the model 
under consideration, the relevant formula is given by 



tot 



01 



rup 



' ^113-^131 



(3.2) 



where Gqi = 7.93 10~^^yr~^ [43] is a precisely calculable phase space factor and e*'^^'^ 
are relative complex phases between the various contributions. These matrix elements 
represent contributions from the direct [M^i^^^^Mxi^^_^\i_^^^) and the neutrino mass {M^) 
mechanisms. 

If mpp is the dominant contribution to OufSfS, using the value of Mi, displayed in 
Appendix 1^, the bound in eq. (|3.1|) can be translated into the limit 



mpp < 460 meV. 



(3.3) 
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If instead the direct contributions are dominant, then assuming the contributions are 
negUgible, eq. leads to [1,2,4,5] 

a;.A'..S^-10-»(^)'. (3.) 

respectively, after taking into account the modifications to the effective Lagrangians dis- 
cussed in Appendix Here AsusY is an effective SUSY breaking scale for the soft terms 
involved in eq. ( |A.3| ), and and Trig/^ are sfermions and gluino/neutralino masses in the 
dominant Feynman diagrams. 

Clearly, the bounds on the RPV couplings depend on the SUSY mass spectra. We thus 
make a simple assumption about the RFC soft SUSY breaking terms: that they follow the 
minimal supergravity (mSUGRA) boundary conditions. The resulting SUSY mass spectra 
are obtained using the spectrum generator SOFTSUSY [16]. At the SUSY scale, all RFV 
couplings are set to zero except for either A'^^^^ or X[^^ and This allows us to emphasise 
effects from the above RFV couplings, without additional complications. Specifically, the 
following set of parameters are defined: 

Mo = [40, 1000] GeV, My2 = [40, 1000] GeV, Aq = 0, tan/3 = 10, sgn// = +1, (3.6) 

where Mq, M1/2 and ^0 are the universal scalar, gaugino, and trilinear soft SUSY breaking 
parameters defined at a unification scale Mx ~ 2.0- lO^^GeV, tan/3 is the ratio of the Higgs 
vacuum expectation values Vu/vd, and sgn^ is the sign of the bi-linear Higgs parameter in 
the superpotential. 
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Figure 6: Ratio TZ = \Mx'^_^^x'^^_^/Mm^i3\, where Mm^p = {rnpp/me)^!^^ in RPV contributions to 
ncutrinolcss double beta decay. The blacked region out at the bottom of the plot is excluded and 
the dotted line delimits regions of different LSP (see text). 

We now turn to interference between different contributions to the Oi^PP rate and 
discuss first the case where both rrif^p and a direct RFV contribution are due to the same 



-7- 



product of RPV couplings. In previous studies, such interference was neglected. For direct 
contributions, this is a good approximation for a SUSY mass scale AgusY of the order 100 
GeV, as the generated from the same RPV couplings is sub-dominant. However, for 
fixed RPV couplings, Af;^/^^ scales as Ag^gY' ^KrsKsi ^'^^^^^ ^SUSY' whereas rripis scales 
as Agygy Thus the rrij^p contribution would dominate at high SUSY breaking scales. The 
ratio 7^ = iMv^^^y^^^M^^^I, where M^^^ = {mpp/me)M^ is shown in fig. ||. The region 
in black at low has either no electroweak symmetry breaking (EWSB), and/or the 

lightest Higgs mass is in violation of the LEP2 direct search limits [44]. The LEP2 95% 
confidence level upper bound implies rrih > 114.4 GeV, but we impose nih > 111.4 GeV in 
order to include a 3 GeV theory uncertainty in the SOFTSUSY prediction of m/j. The yellow 
dotted line on the left separates the region with a stau (left) and a neutralino (right) LSP. 
We see that while in the lower mass region the direct contributions dominate over the m^^ 
contributions, they become comparable in the high mass region, where interference cannot 
be neglected. We thus include interference terms in the calculations which follow. 

On the other hand, despite the Ag^gy dependence on My^^^, the ratio |Mv^^^/Mm^^| 
is greater than 20 in our parameter space region. This is because generated by X'm 
is heavily suppressed by the running up quark mass insertions in the loop diagrams. For 
this reason, we may safely neglect m^^ contributions to Ovfip when considering cases with 
non-zero X'm unless it originates from some other coupling. 

4. B^-B^ mixing, single slepton production at the LHC and 0^/3(3 
4.1 Implications of A'^]^3A']^3]^ bound from B^-B^ on Oi//3/3 

It was shown in [4] that for sparticle masses of ^ 100 GeV, a stringent bound on A'^^^gA'^^g^^ 
comes from O^pp. Another competing bound comes from B^-B^ mixing. A Feynman 
diagram of the latter process is displayed in fig. 0. A recent update [45] by the CKMfitter 




Figure 7: B^-B^ mixing through coupling product X'njX'i^i 



group shows that at 95% confidence level, the magnitude of any new physics effect to B^-B^ 
mixing must be less than the SM contribution. This result is shown in fig. In this figure, 
the Standard Model solution is located at = 1, and deviation from unity represents new 
physics effects contributing to B^-B^ mixing. The upper limit of A'^^^gA'^gj^ is obtained by 
updating the results in [4] to take into account the latest B^-B^ mixing data. We obtain 
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Figure 8: Possible new physics contributions to B'^-B^, from Ref. [66]. The 95% C.L. region is 
shown in soUd orange Hne. The SM solution is located at A^^ = 1, and deviation from the SM value 
may be attributed to the RPV contributions proportional to A'^j3A'j3j. 



To compare limits on X'n^X'i^^ from B^-B^ and from {)v(3f3, we recall that the bound 
from Oiypp depends on the sbottom mass squared matrix. In the case where all SUSY 
breaking parameters are of the same order of magnitude, this bound relaxes approximately 
as the cube of the sbottom mass scale, which is more rapid compared with the B^-B^ bound 
in eq. ([4.1|). However in the most general MSSM the mass parameters relevant for these 
two bounds are independent. Which of the bounds is more stringent depends therefore on 
the ratio of sneutrino to sbottom masses. 

In the following, we restrict our discussion to the parameter space discussed in section 
|3[ From [4], we expect the Oz//3/9 bound to be more stringent only in very low AsusY of 
around 100 GeV. With the new B^-B^ limit, and the fact that mSUGRA-like mass spectra 
generally have squarks much heavier than the sleptons, we find that the bound from B^-B^ 
mixing is more stringent than that from OvP/S in all allowed regions of parameter space we 
explore. This means that if the direct X'n^X'i^i contribution is the only source of Oz^/?/?, then 
Ty^^ must be larger than the current experimental limit. This conclusion is sensitive to 

theoretical uncertainties in our predicted value of Ty2^ coming from the NMEs. Estimates 
in such uncertainties vary: NME calculations based on different nuclear model assumptions 
and input parameters could differ by a factor of 1.3, 3 or even up to 5 [64,65]. These are 
then squared in order to obtain the half-life prediction. We take for example a factor of 
3 uncertainty in the predicted half life, equivalent to ±0.5 in logiQT^J^^ . If the NMEs 

predicted Ty^^ to be a factor of 3 less than the values taken here, then there is a small 

region where the stronger bound would instead originate from T^^^ at the lowest viable 
values of M1/2 and Mo/GeV ~ 400 - 600. See the appendix for details about the NMEs. 

The variation of the T^^^ lower limit, obtained from an upper bound on A'^^^gA'^^g^^ from 
B^-B^ mixing is shown in fig. |9|. As in fig. the (blue) dotted line separates regions with 
stau LSP and neutralino LSP, while the black region at the bottom is excluded due to no 
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Figure 9: Lower limit on T'[1l^'\'^^Ge), using upper limit on X'n'i^'i^i obtained from B^-B^ mixing. 

The black contour shows where a T^j2^ limit of 10^^ years is expected. The black region at the 
bottom of the plot is excluded and the dotted line delimits regions of different LSP (see text). 

EWSB or the higgs being too light. We see that Ty^'^ ~ lO^'^-lO^^yrs is still allowed in 
much of the parameter space, so that Oi^/5/3 can be detected by next generation experiments. 
In particular, there exist good prospects of observing a Oiyf3f3 signal in the region with 
relatively low Mi^2- This is because the sbottom masses receive large renormalisation 
effects from the gluino mass and as a result are much lighter in the low region. This 
in turn enhances ^A'^^gA'^g^ compared with corresponding values in the high M1/2 region. 




1.9e+25 le+26 le+27 1.9e+25 le+26 le+27 

observed Ge half life (yrs) observed Ge half Ufe (yrs) 



Figure 10: Effect of a near-future measurement of T^J^'^ for Mq = 680 GeV, M1/2 = 440 GeV, 
given current B^-B^ mixing constraints. In both panels, the shaded regions are allowed. 

We now consider the case where sbottom exchange plus other possible processes con- 
tribute to Ou(3(3, and where m^^ may not be solely due to y^i^i^n^ 7^ leading to the 
process m fig. |. Thus, we consider that there may be other contributions to m^^, coming 
from bi-linear RPV couplings, or Ajj^ couplings, for example. We imagine that LHC mea- 
surements are compatible with lepton-number violating mSUGRA signals, with Mq = 680 
GeV and M1/2 = 440 GeV, ^0 = and tan/3 = 10. In practice, these numbers would 
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be determined with some uncertainties, which we ignore for now since we are merely illus 

1/ 



trating the point. Then, a measurement of T^^^^ in the next generation of experiments 



could constrain the Qvj3j3 mechanism when combined with B^-B^ constraints. In order to 
quantify how much of the Ou(5(5 width may come from direct processes involving sbottoms, 
we define 



My A' 



Mi 



tot 



(4.2) 



where Mtot is the total matrix element including both sbottom mediated and m^^-induced 
contributions, r = 1 implies that the sbottom-mediated contributions could account for 
all of the Qu(5(5, whereas r < 1 requires some extra component (for example from m^p). 
r > 1 also requires an additional component that destructively interferes with the sbottom 



exchange process. We illustrate this in fig. 10, where the hatched region in the left-hand 



panel shows which values of r would be allowed by current B^-B^ mixing constraints. 
We see that, for T^"!^^ < 10^^ years, it is not possible to explain with sbottom 

exchange while simultaneously satisfying B^-B^ bounds. The right-hand panel shows what 
the inferred value of mpp may be assuming that the only contributions are from sbottom 
exchange and m^^. The range of values comes from the possible size of the direct sbottom- 
mediated contributions and the fact that the interference could be either constructive or 
destructive. We again see that T^^^ < 10^^ years would imply that the direct contribution 

may not account for Ou^P alone. Theoretical uncertainties in the T^^!^^ predictions coming 
from NMEs would affect the inferences in fig. widening the band in the right-hand panel 
and raising the bound in the left-hand panel. 

4.2 Single selectron production at the LHC via A'^^^^^ 




Figure 11: Production of a single selectron at resonance via X'm, followed by a gauge decay into 
an electron with a neutralino LSP. The LSP further decays into three final states via a virtual 
sparticle, leading to a same sign, di-lepton signal for the whole resonance production process. 



Depending on the value of A'^^^^^, resonant production of a single selectron may be 
observed at the LHC. The related process of single smuon production is studied in [46] where 
like sign di-lepton signals are used because of small backgrounds. A diagram showing like 



sign di-lepton production, with decay of a selectron via a neutralino is displayed in fig. 11 



Earlier studies based on different signatures can be found in [47-51]. For previous analyses 
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of lepton number violation utilizing the same sign di-lepton signature see e.g. [52-55]. A 
first study of single slepton production in stau LSP scenarios can be found in [56]. However, 
to the best of our knowledge, the discovery reach of such stau LSP scenarios at the LHC is 
not available in the literature, and hence our following discussion will be restricted to the 
case with a neutralino LSP. 

At low sparticle mass scales of ~ 100 GeV, the stringent bound from Quf3(3 renders 
single slepton production unobservable. However the strong dependence of this bound on 
the SUSY mass scale, as shown in eq. ( |3.5| ), means that it may be possible for this process 
to be observed at higher SUSY mass scales. 




200 400 600 800 1000 



Mo/GeV 

Figure 12: mSUGRA parameter space in which single slepton production may be observed at the 
LHC for mpfi = 0, tan/3 = 10, = and lOfb^^ of integrated luminosity at 14 TeV centre of 
mass energy. In the top left-hand black triangle, the stau is the LSP, a case not covered by this 
analysis. The bottom black region is ruled out by direct search constraints. The labelled contours 
are extracted from Ref. [46], and show the search reach given by the labelled value of X'm- From 
bottom to top, the white, dark-shaded and light-shaded regions show that observation of single 
slepton production at the 5cr level would imply T°J^^ < 1.9 • lO^^yrs, 100 > T°J^^ /lO^^yvs > 1.9 
and T^j2^ > 1 x lO^^yrs, respectively. The upper and lower dashed curves show where the contour 
between the dark-shaded and light-shaded regions would move to if [m/j/jl = 0.05 eV were included 
with constructive or destructive interference, respectively. 

We will now compare the discovery reach of A211 [46] at the LHC via single smuon 
production with the bounds on X'^^^^ coming from Oz//3/3. As Ref. [46] does not include 
detector effects, the reach of X'211 is expected to be readily applicable to X'^n without 
significant changes. For this reason we shall use the discovery reach of A211 interchangeably 
with that of A'^^^^ from now on and use the results of Ref. [46] as an estimate of the 5a- 
significance discovery reach for 10 fb^^ of LHC integrated luminosity at a centre of mass 
energy of 14 TeV. Cuts were placed on the leptons: a minimum transverse energy cut and 
an isolation cut in order to reduce heavy quark backgrounds. Other SUSY processes are 
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cut by requiring a maximum missing transverse energy and requiring that there are at 
most 2 or 3 jets above a minimum transverse momentum. We refer the interested reader 
to Ref. [46] for more details. 

Fig. |l2| shows regions of the Mq — M1/2 plane where single slepton production may be 
observed via like-sign electrons plus two jets. The black regions with high lo'^ 
Mq have a stau LSP. As discussed before, this is not included in our discussion because 
a detailed Monte Carlo study is not yet available (however see [56] for an initial study). 
The black regions with small excluded as in fig. |3|. In the white region, single 

slepton production by X'm could not be observed without violating the current bound 
upon Ty^^ . The darker shaded region shows where the observation of single slepton 
production at 5a above background implies that {)v(3f3 is within the reach of the next 
generation of experiments, which should be able to probe T^^^ < 1 x lO^^yrs [57,58]. 
Conversely, if Oz//?/3 is discovered by the next generation of experiments, we should expect 
single slepton production to be observable and test the X'm hypothesis. We do not expect 
Aq or tan P to affect the shape of the regions much, since they have a negligible effect on the 
selectron mass and the couplings in the relevant Feynman diagrams. In the light shaded 
(upper) region, a 5a single slepton discovery at the LHC implies that the next generation 
of experiments would not be able to observe Qv(5(5. Conversely, if Qv(3(3 is within reach of 
the next generation of experiments, the LHC would see single slepton production signal in 
this region at greater than 5a significance. 

If Qvf3f3 is marginally observed at the limit of the next round of Qv(5(3 experiments, 
it is possible that a contribution from mpp could contribute significantly (if it originated 
from a different coupling to A'^^^^). For ''^Ge, a reach of Ty!^^ ~ 10^''yrs implies that an 
inverted or quasi-degenerate neutrino mass spectrum may contribute to at an ob- 

servable level. Whether it interferes with the direct contribution amplitude constructively 
or not will affect the potential observability at the LHC, and should be taken into account. 
We shall assume that |m/3/3| = yAm^^ ~ 0.05 eV, as implied by neutrino oscillation 
data. We show in fig. 12 where the upper edge of the darker region would move to and 
constructive (upper dashed curve) or destructive (lower dashed curve) interference between 
the mpj3 contribution and the direct contributions. We see that for the case of constructive 
interference in fig. [l^, discovery of single slepton production becomes very difficult with 
10fb~^ of data. On the other hand, the situation improves substantially if destructive 
interference occurs instead. To see this, we note that for a fixed t'^J^^ , the introduction 
of a non- vanishing m^p which interferes destructively with My^^^ would imply an increase 
in the direct contributions. This would require an increase in A'^^^^, or a decrease in the 
mass of the SUSY spectrum, or both, all of which lead to an increase in single slepton 
production rate. We see that single selectron production may be observed for ^ 500 
GeV. If the neutrino mass spectrum is instead normal hierarchical, the diagram will 
be sub-dominant for the half life discussed above. This situation may be approximated by 
setting mpis = 0. 

We show in fig. 13 the variation of the discovery reach of A'^^^^^ with Mq along the 
line Ml/2 = 300 GeV + O.6M0. Above the dotted light line, single slepton production 
will be observed at the LHC. We see from the figure that for nearly all of the parameter 
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Figure 13: Comparison o{T^j2^ and single slepton discovery reach as a function of A'j^j^j^ along the 
mSUGRA slope Afi/2 = 300 GeV + 0.6Mo, with = 0, tan/3 = 10 and sgn^ = 1. The black region 
on the top left corner is ruled out by 0iy(3/3. The region above the solid black line is accessible in 
near future Qv^P experiments, whereas the light dotted line shows the lower limit of X'm for single 
slepton production to be discoverable at the LHC. 



space where Ou(3(3 can be measured by the next generation of experiments, the LHC would 
provide a confirmation of the supersymmetric origin of the signal by observing single slepton 
production at the 5a level. 
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Figure 14: Effect of a near-future measurement of T'^J^'^ for Mq = 680 GeV, M1/2 = 440 GeV, 
given a 5a observation of single slepton production at the LHC in 10 fb^^ at 14 TeV centre of mass 
energy. In the right-hand panel, the shaded region would be allowed. 



We now ask the question: if single slepton production were observed at the LHC, and 
a measurement of T^^^^ were made, what could be divined about the relative contributions 
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between direct or Majorana neutrino-induced Oi^/3/3? We define 

Mx' 



r' 



Mi 



tot 



(4.3) 



where Mtot is the total matrix element including both sbottom mediated and m^^-induced 
contributions. At our parameter test point of Mq = 680 GeV, M1/2 = 440 GeV, ^0 = 
and tan /3 = 10, a measurement of the single slepton production cross-section could be 
used to infer a value of X'ln, which would then imply a value of r'. Here, we assume that 
the cross-section was just at the 5a level above background and show which r' could be 



inferred from the measured value of T^!^^ in the left-hand panel of fig. 14. The figure shows 



a definite prediction for r', which is 1 at a half-life of 10^^ years, where Qvj3j3 could come 
entirely from the direct contribution. In general, r' > 1 requires destructive interference 
between the X'm contribution and another contribution (either A'^3]^A'^]^3 or nipp). Assuming 
the relative phases between different contributions are real, in the destructive interference 
case there are two solutions to I-^a'^^^ ~ -^otherl ~ l-^iot|j where \Mtot\ is a prescribed 
constant derived from Ty!^^ and eq. (^.2|). We note that NME uncertainties, which are 
neglected here, would turn the definite prediction into a band of possible predictions. The 
r' prediction would acquire error bands, at the level of ~ 0(20)%, from measurement errors 
in the SUSY spectrum, cross-section and NME uncertainties. A further dedicated study 
including simulations of the experiments is required in order to quantify these latter errors 
more exactly. 



The right-hand panel of fig. 14 shows correspondingly what may be deduced about 
m^^, also neglecting measurement and NME errors, but taking into account two possible 
contributions to the matrix element: My and Mx'_^^^- If instead we assumed that {)v[3[3 
was due entirely to the Majorana neutrino induced contribution, then one deduces m^p 
from Ty^^ as in fig. |5|. We see that the definite prediction widens to a band, due to the 
relative complex phase between the two contributions to the matrix element. Any inferred 
value m^^ > 10 meV corresponds to a non-hierarchical pattern (i.e. inverted or quasi- 
degenerate) of neutrino masses, and so one would infer that the hierarchy is non-hierarchical 
for Ty^^ > 10^'' years. We see that, for this parameter space point, if Ty^^ 10^^ 
years, the observation of single slepton would imply that a normal hierarchy is still viable, 
although still one could not tell without additional measurements if the neutrino mass 
hierarchy were inverted or hierarchical. On the other hand, for Ty!^^ > 1.1 x 10^^ years 

or T^j2^ < 0.9 X 10^^ years, neutrino masses have the non-hierarchical pattern. It remains 
to be seen how large this window remains after measurement uncertainties are calculated. 
In the hmit where the contributions from m^^ and My are real, a given r°/f ^ would 
result in two-fold possible predictions: either the top or the bottom of the bands in the 
right-hand panel of fig. 0. Theoretical errors originating from the NME calcuations would 
widen the band further. 

5. Discussion and summary 

In this paper, we have discussed the interplay between a number of observables in RPV 



- 15 - 



SUSY: Oiy(3(3, B^-B^ mixing and first generation single slepton production at the LHC. 
We saw that, while it is difficult to infer unambiguously the presence of RPV effects from 
individual observables, further insight could be gained if these information are analysed 
simultaneously. Bounds from B^-B^ mixing and single slepton production could constrain 
the extent to which X'n^^i^i and A'^^^A'^^i operators may contribute to Quf3f3. They provide 
extra handles to complement other strategies for divining the underlying Qvf3f3 mechanisms, 
for example by comparing life times of different nuclei [27]. 

There are a couple of caveats one should bear in mind however. Firstly, our numerical 
inference between collider observables and Oi^/?/? depends on the relation between the masses 
of the SUSY particles that mediate the different processes. Constraints on these masses 
from experiment are required: ideally directly, but otherwise because some other simple 
theoretical SUSY breaking ansatz such as mSUGRA fits LHC data well. Secondly, the 
estimated T^!^^ is subject to uncertainties in the NME. This is important particularly 
when using B^-B^ data because the sensitivity of the next round of '^^Gc Oi^/3/3 experiments 
is not expected to go far beyond the 10^^ yrs level. The bounds from Oz//3/3 and B^-B^ are 
comparable in regions with 0.1-1 TeV scale SUSY breaking, with the latter being slightly 
more stringent in the parameter space we explored. An increase in the NME value used 
might change this observation. 

At present, B^-B^ mixing implies only an upper bound on the possible contribution 
of heavy sbottom exchange to any measurement of Ty^^ , given information on the SUSY 

spectrum from the LHC. Sbottom exchange could be solely responsible for T^J^^ (^^Ge) 
of ~ lO^^yrs, which is potentially observable in the near future. However, the measured 
value of Ty^^ could still require there to be a non-zero Majorana neutrino contribution, 
depending upon its value. A much more informative two-sided bound on the sbottom 
exchange contribution would result if a future measurement of B^-B^ mixing required an 
extra component from outside the Standard Model. 

We point out that, contrary to previous expectations, a same sign di-lepton signal 
from single selectron production via gauge decay at the LHC could be observed. This is 
despite the stringent bound from Oz//3/3, because the constraint from Ov^P on A'^j^;^ relaxes 
rapidly as the SUSY scale increases. We have shown that if a direct contribution via A^^^^ 
is the dominant Oi^/3/3 mechanism and OuP/d of ^^Ge is just beyond the current reach, there 
is a good chance of observing single selectron production at the LHC. Such a scenario is 
not ruled out by current Qvf3f3 bounds for a heavy enough SUSY spectrum. Knowledge 
of the SUSY spectrum can be combined with OvfSf] and single-selectron data to bound 
the X'li^ contribution to Ov^p. Thus, evidence for other contributions may be obtained 
(for example from Majorana neutrino exchange) and the size of the other contributions 
bounded. Under the hypothesis that only the A'^^^^ and nifj/j contributions are significant, 
mpp may be deduced and information about the neutrino mass spectrum is thus obtained. 
For some ranges of T^^^^ , this could settle the question of whether the neutrino spectrum 
were hierarchical or not withoTit the inclusion of other observables. 

It will be an interesting exercise in future work to examine a particular point in param- 
eter space with a dedicated LHC simulation study in order to quantify the errors obtained 
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on the inferred contributions to Qv(3(3 in a combined fit. 
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A. Parton Level Contributions to Oi^/3/3 

A.l Light Majorana neutrino exchange: Qv(5(3 via rapp 

The effective Lagrangian after integrating out the W gauge boson, and the ALe = 2 
Lagrangian with a virtual Majorana neutrino are 

Gf 



C 



'-EW\^I 

eff,AL,=2, 
EW 



V2 



67^(1 - 75)z^ UyJ^{l - 



67^(1 -75)^-7^^6 



(A.l) 



respectively, y is a colour index, and in the second line of eq. ( [A.l ), Jy_^ = tiy7'^(l — 75)^^. 
A Feynman diagram depicting this process is displayed in fig. ^. 

A. 2 Heavy sbottom exchange: via X'u^X'i^i 

In the basis where both the down-type quark and the charged lepton mass matrices are 
diagonal, the coupling product A'^]^3A']^3]^ leads to an effective Lagrangian involving exchange 
of one SUSY particle of the form 



-^113-^131 



^v'^iUhiNsU + 75)e'){uya,4l + 75)^^) 

+2(i?i(l + 75)e'=)(7i,(l+75)d^ 



(A.2) 



where [59]^ 



E 



KmkKkli^CKM)lm 
2V2GF 



m- m- 
"fc(i) 



■fe(2) 



-^n3-^ra31 _ 

OLL ban 



i>LR 



(A.3) 



^The first term in eq. (A.2) diflers from [4,5] and aiso tfie preprint version of [59] by a factor of i. Our 
check agrees with the published version of the latter reference. 
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In eq. ( [A.2| ), Upmns is the 3x3 unitary PMNS neutrino mixing matrix [60,61], and 



a^^^ = ^7^,7'^]. In eq. (|A.3D, Uckm is the CKM matrix and m?- , m?- and m?- 

<ik{LL) d.k{LR) d.^(RR) 



denote entries in the k-th generation down type squark mass squared matrix. In particular, 
m^^ and m^^ denote the 2 sbottom mass eigenvalues, and is the sbottom left-right mixing 
angle. The relations between the mixing angle and the entries in the mass and flavour basis 
sbottom mass matrices follow those in SOFTSUSY [16]. 

The complete Ovf^/S matrix element with a leptonic current coupled to a quark current 
via a W boson does not contain a neutrino mass insertion, and hence is not suppressed by 
the light neutrino mass scale. The ALg = 2 Lagrangian is given by 

£^'v;;fSi.,w = f v[i(.->(i-7.)ie=)jr,jp. 

where Jps = Uy(l + 75)^^ and J^*^ = Uya^'^{l -|- 75)^^ are the pseudo-scalar and tensor 
currents respectively. A Feynman diagram depicting this process is displayed in fig. ^. The 
matrix element is [4, 5] 

My^^^y^^^=rj\Mf + M^), (A.5) 

rj^ is defined as in eq. ( |A.3| ) with n = 1, and and AI"^ denote the 2 nucleon mode and 
pion mode contributions and will be detailed in Appendix 

A. 3 Sparticle exchange: via X'mX[i^ 

Following the notation of [1], the effective Lagrangian with X'^^^ in the direct RPV Gvf3[3 
process involving exchange of two SUSY particles is given by 

^I(f;Al'r'(-) = ^-;ne-(i+75)e^] 



X 



(^9 + Vx)(JPsJpS - ^Jt'^Jt^lu) + {Vxe + ri'g + T] ?)JpsJps\ (A.6) 



where the RPV coefficients are defined to be 



vros X'lii rUp 1 1 



2 



6 rria m% m- lm:~ m 

^ 7Ta2 Xf^^ rUp e%{u) e%{d) eLAu)^RM) ■ 
"Ix 2 Gl ^ m.^ ^ m% ^ m% ml w?- 



Gi mv. mt 



F j=i Xi ei 



, _ 27ras mp 1 
~ o ^2 V 



Afii A rUp eijujepj^ gL,MgL,(e) eije)epjd) 
7^? = vra2^y5- / V 9 9 2 9 2 2 J' l^-'J 

F j=l Xt UL dR EL eL dR 
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and again we follow the notation of [1]. The e's denote rotations between mass and gauge 
eigenbasis in the gaugino-fermion-sfermion vertices. To facilitate comparisons with the 
literature, we display only the first generation sparticle contribution above but include 
contributions from all three generations in the numerical calculations. The relevant Feyn- 
man diagrams are displayed in fig. ^. Note that our expressions above are different from 
those presented in [1], [2] and [62]. Reference [1] made an approximation when extracting 
the colour singlet currents from the Feynman diagrams. Also, the expressions for rj^ and 
rj^^ (following the convention of [1]) are different, with the discrepancy coming from the 
colour flows of fig. 4(e) and fig. [4(f)| . References [2, 62] took proper account of the colour 
singlet extraction, which we have checked independently. However rj^ and r]^j remained 
the same as in [1]. In the parameter space we explore, the numerical differences induced 
by such changes in the r]s are small compared to the uncertainty in the nuclear matrix 
elements. 

Mx[,, is given by [1,2] 

Mx' = (rig + r/;,)M|^ + (t?^, + r]', + %f)Mf 



+ 8 [{v-g + Vx) + ^{Vg + Vx+ Vxe + v'g + V^f)) {^M^^ + M^^). (A.8) 

The rys correspond to those defined in eq. (A^), and -^|^) ^^'^ ^'^^ M^'^ denote the 2 
nucleon, 1 pion and 2 pion exchange modes respectively. They are detailed in Appendix^. 



B. Nuclear Matrix Elements 



channel 


NME(^6Ge) 


value 


Ref. 


^'iii-^ni 




283 


[1] 




Mf 


13.2 


[1] 






-18.2 


[2] 






-601 


[2] 


'^113'^131 




-0.9 


[5] 






604 


[5] 






2.8 


[63] 



Table 1: Nuclear matrix elements of Ge used. The value of the 2 nucleon mode contribution to 
M\' X' includes a factor of i discussed in footnoted. For model details of the NME calculations, 
we refer readers to the literature. 



As we are primarily interested in collider effects due to couplings from the direct 
contributions, we adopt numerical values collected from various sources in the literature 
for our discussions here. These parameters are displayed in table |l|. 
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